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1. Introduction 

Compton scattering at low energies is an invaluable tool with which to study the electromag- 
netic structure of hadrons. The real Compton scattering amplitude describing the elastic scattering 
of a photon on a spin-half target such as the proton or neutron can be parameterised as 

TyN = Ait' ■s+A 2 e' -kt-k! + z'A 3 (7 • (e' x e) + iA A o -(k 1 xk)t' t 

+iA 5 d- [(t' xk)t-k'-(txk')t' -k] +iA 6 d- [(f x k')t-k'- (e x k)t' ■%] , (1.1) 

where we work in the Breit frame of the system and the incoming and outgoing photons have 
momenta k = (co,k = cok) and k' = (co,k' = cok'), and polarisation vectors e and e' , respectively. 
The A,- = A, (ft), 0) are scalar functions of the photon energy and scattering angle, cos = k - k , 
and can be separated into two pieces. The first set, the Born terms, describe the interaction of the 
photon with a point-like target with mass, M^, charge, eZ (where e > 0), and magnetic moment, /x. 
The remaining parts of the amplitude describe the structural response of the target and expanding 
the amplitude for small photon energies relative to the target mass and keeping terms to &{(o i ) one 
can describe the target structure in terms of the electric, magnetic and four spin polarisabilities [0], 
a, j6, and 71-4, respectively. 

Lattice techniques provide a method to investigate the non-perturbative structure of hadrons 
directly from QCD. In particular, the various hadron polarisabilities can be computed. Direct calcu- 
lations of the required hadronic current-current correlators are difficult on the lattice and so far have 
not been attempted. However progress has been made ^J,|5|] in extracting the electric and mag- 
netic polarisabilities by studying the quadratic shift in the hadron mass that is induced in quenched 
lattice calculations in constant background electric or magnetic fields. These studies have inves- 
tigated the electric polarisabilities of various neutral hadrons (in particular, the uncharged vector 
mesons and uncharged octet and decuplet baryons), and the magnetic polarisabilities of the baryon 
octet and decuplet, as well as those of the non-singlet pseudo-scalar and vector mesons. As we 
shall discuss below, generalisations of these methods using non-constant fields allow the extraction 
of the spin polarisabilities from spin-dependent correlation functions and also allow the electric 
polarisabilities to be determined for charged hadrons. More generally, higher-order polarisabilities 
and generalised polarisabilities are accessible using this technique. 

As with all current lattice results, these calculations have a number of limitations and so are not 
physical predictions that can be directly compared to experiment. For the foreseeable future, lattice 
QCD calculations will necessarily use quark masses that are larger than those in nature because of 
limitations in the available algorithms and computational power. Additionally, the volumes and 
lattice-spacings used in these calculations will always be finite and non-vanishing, respectively. 
For sufficiently small masses and large volumes, the effects of these approximations can be inves- 
tigated systematically using the effective field theory of the low energy dynamics of QCD, chiral 
perturbation theory (#PT) In this paper we present results for the nucleon electromagnetic and 
spin polarisabilities at next-to-leading order (NLO) in the chiral expansion. We do so to discuss 
the infrared effects of the quark masses and finite volume in two-flavour QCD and its quenched 
and partially-quenched analogues (QQCD and PQQCD). The polarisabilities are particularly inter- 
esting in this regard since they are very sensitive to infrared physics and their mass and volume 
dependence is considerably stronger than that expected for hadron masses and magnetic moments. 
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2. Polarisabilities on the lattice 

The use of external fields in lattice QCD has a long history. The pioneering calculations of 
Refs. [0, |, g [□]] attempted to measure the nucleon axial couplings, magnetic moments and 



electric dipole moments by measuring the linear shift in the hadron energy as a function of an 
applied external weak or electromagnetic field. As discussed in the Introduction, various groups 
[Q, |5|] have also used this approach to extract electric and magnetic polarisabilities in quenched 
QCD by measuring a quadratic shift in the hadron energy in external electric and magnetic fields. 
The method is not limited to electroweak external fields and can be used to extract many matrix 
elements such as those that determine the moments of parton distributions and the total quark 



contribution to the spin of the proton [|12j]. Here we focus on the electromagnetic case. 

The Euclidean space (X4 = x) effective action, 5 e ff[A] = / d 3 xdx J5f e ff (x, T;A), describing the 
gauge and parity invariant interactions of a non-relativistic spin-half hadron of mass M and charge 
q with a classical U(l) gauge field, A^(x, x), is formed from the Lagrangian 



J^ eff (x,T;A) = ^(JC.T 



-^ + iqA 4 ^j + <y lV m qA ^ -pLC-H + 27i(aE 2 -fiH 2 ^ (2.1) 



-2th" (-7£ l£l a • E x E + y Ml M, o-HxH + y MlE2 a i E ij H i + y ElM2 o>H ij E 



¥(x,t) + . 



where E = — 4-A(x, x) — VA^(x, x) and H = V X A(x, x) are the corresponding electric and mag- 
netic fields, X = -^X denotes the Euclidean time derivative, X'-i = ^(d'X J : + d->X'), and the ellip- 



sis denotes terms involving higher dimensional operators. The constants that appear in Eq. ( |2.i| ) 
are the magnetic moment and electromagnetic and multipole polarisabilities: Ye^Ei = — (/i +73). 
YMiMt = 74, Ye x m 2 = 73 and YMiE 2 = 72 + 74- The Schrodinger equation corresponding to Eq. (jD]) 
determines the energy of the particle in an external U(l) field in terms of the charge, magnetic 
moment, and polarisabilities. 

Lattice calculations of the energy of a hadron in an external U(l) field are straight-forward. 
One measures the behaviour of the usual two-point correlator on an ensemble of gauge configura- 
tions generated in the presence of the external field. This changes the Boltzmann weight used in 
selecting the field configurations from det [p+ m] exp [S g ] to det \p+ i Qfi + m] exp [S g ] , where 
p is the SU(3) gauge covariant derivative, Q is the quark electromagnetic charge operator, and S g is 
the usual SU(3) gauge action. Since calculations are required at a number of different values of the 
field strength in order to correctly identify shifts in energy from the external field, this is a relatively 
demanding computational task. The exploratory studies of Refs. @ [| ||] used quenched QCD in 
which the gluon configurations do not feel the presence of the U(l) field as the quark determinant 
is absent. In this case, the external field can be applied after the gauge configurations had been gen- 
erated and is simply implemented by multiplying the SU(3) gauge links of each configuration by 
link variables corresponding to the fixed external field: (x)} — > (x)exp[ieaA^]}, where 
a is the lattice spacing. These studies are interesting in that they provide a proof of the method, but 
the values of the polarisabilities extracted have no connection to those measured in experiment. 

It is clear from Eq. ( |Q| ) that all six polarisabilities can be extracted using suitable space and 
time varying background fields if the shift of the hadron energy at second order in the strength of the 
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field can be determined. In order to determine the polarisabilities, we consider lattice calculations 
of the two-point correlation function 

C ss ,(p, x;A) = J d 3 xe^(0\x s (x, t)^(0,0)|0) a , (2.2) 

where % s {x, x) is an interpolating field with the quantum numbers of the hadron under consideration 
(we will focus on the nucleons) with z component of spin, s, and the correlator is evaluated on the 
ensemble of gauge configurations generated with the external field, A' 1 . 



For weak external fields (such that higher order terms in Eq. (2.1) can be safely neglected), the 



small p and large x dependence of this QCD correlation function is reproduced by the equivalent 



correlator calculated in the effective theory corresponding to the Lagrangian, Eq. (2.1). That is 

C ss ,(p,x;A) = J d 3 xe'^^j-r J SH* W^(x,T)^,(0,0)exp (-S eff [A]) , (2.3) 
where Jeff [A] = / W f Wexp (-S e ff [A]). Since tne right-hand side of Eq. (O) is completely 



determined in terms of the charge, magnetic moment and polarisabilities that we seek to extract, 
fitting lattice calculations of C ss /(p,x;A) in a given external field to the effective field theory ex- 
pression will enable us to determine the appropriate polarisabilities. In the above equation we have 
assumed that the ground state hadron dominates the correlator at the relevant times. For weak fields 
this will be the case. However one can consider additional terms in the effective Lagrangian that 
describe the low excitations of the hadron spectrum. 

In many simple cases such as constant or plane-wave external fields, the EFT version of 



C ss '(p,x;A) can be determined analytically in the infinite volume, continuum limit [14]. How- 
ever at finite lattice spacing and at finite volume, calculating C ss i(p, x;A) in the EFT becomes more 
complicated. In order to determine the EFT correlator, we must invert the matrix Jff defined by 

5ia tt [A] = £Y,L w U^xWsAx,ix,y,v>MVs(y,^) , (2.4) 

where Si att [A] is a discretisation of the EFT action in which derivatives are replaced by finite differ- 
ences. For the most general space-time varying external field, this must be inverted numerically; 



given a set of lattice results for the correlator, Eq. (2.3) is repeatedly evaluated for varying values 
of the polarisabilities until a good description of the lattice data is obtained. If the external fields 
are weak, |A^(x, x)\ 2 <C Aq CD for all x and X, a perturbative expansion of J^ _1 in powers of the 
field can also be used. 

To extract all six polarisabilities using such an analysis, we need to consider a number of 
different fields; lattice calculations of the correlators in Eq. (also for different spin config- 
urations) using the exemplar fields given in Ref. ^ for various field strengths are sufficient to 
determine the full set of polarisabilities. As an example, the behaviour of the correlator in the field 
A^j(x) = (iaiX, 0,0,0) (which corresponds to a constant electric field in the xi direction) is given 
by 

C s/ (p,x;A w ) ~ 5 v exp{-^| [ ai {q 2 x 2 + \2M%a) -3iqX Pl ] } e - Mx e -^ + 0{t%) 



— > <V ex P 



-(M + 2W)t-^1t 3 
v 11 6M 



+ (9 {a}) . (2.5) 
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In this case, the perturbative series has been resummed exactly in the continuum, infinite volume 
limit and the higher order corrections come from terms omitted in Eq. (2T ). For electrically neutral 
particles, the exponential fall-off of this correlator determines the polarisability a once the mass 
M has been measured in the zero-field case. When a charged particle is placed in such a field 
it undergoes continuous acceleration in the x\ direction (this is described by the T 3 term in the 
exponent). However at times small compared to \f6M/qa\, the correlator essentially falls off 
exponentially. Matching the behaviour of Eq. ( p3| ) to lattice data for a charged hadron will again 
enable us to determine the electric polarisability, a. Whilst the above analysis assumed infinite 
extent in the x\ direction, it remains valid for L <C m^ 1 . 



3. Finite volume effects in nucleon polarisabilities 

To calculate the quark mass and volume dependence of the nucleon polarisabilities, we use 
heavy baryon chiral perturbation theory. Our results also include the quenched, and partially- 
quenched versions of this theory appropriate for existing lattice calculations. The details of the 
calculations are given in Ref. [|T|], here we merely summarise the result for the case of the fourth 
spin polarisability in the QCD case. It is convenient to separate the different contributions as 



74 



7T maly + tf + 7 1 4 00p , (3-D 



corresponding to the contributions from the anomalous decay of flavour neutral mesons to two 
photons, Born-terms involving the A-isobar resonance, and loop diagrams, respectively. 

The anomalous decay of flavour neutral mesons to two photons has important consequences 
in Compton scattering in non-forward directions. Anomalous decays are well understood in %?T, 



entering through the Wess-Zumino-Witten Lagrangian [15J; extensions to the partially-quenched 
case are discussed in Ref. |jl|]. The contributions to the amplitude from the Born-terms involving 
the A-isobar resonance are also significant. The resultant terms are found to be 

anomaly 3e Ganom A & /o t\ 

74 167T 3 / 2 m 2 ' 74 ~ 2167r(2M;v) 2 A 2 ' 1 ' 

where G anom = £a(2Z — 1) and/^ is the magnetic dipole transition coupling |JIJ]. 

For the loop contribution, using the effective couplings Gb = \g\ and Gj = qgj^ (which are 
modified in the quenched and partially-quenched theories [[T]]), we find that [0: 



i e 2 

>op ^ e 



1927T 3 / 2 



G B G T 1 A 

+ -^< T ~t~ TTTn tTTTT 



4m 2 3 |A 2 -m 2 2(A 2 -m 2 ) 3 / 2 



A— a/A 2 — m 2 + ie 



A + a/A 2 — m 2 + ie 



.(3.3) 



The finite volume of a lattice simulation restricts the available momentum modes and conse- 
quently the results differ from their infinite volume values. Here we shall consider a hyper-cubic 
box of dimensions L 3 xT with T ^>L and m n L 3> 1 (smaller volumes in which m n L ~ 1 [lM 19] are 



also discussed in Ref. [jl]]) which leads to quantised momenta k = (ko,k), k = = ^(ji,j2)73) 
with jj G Z, but ko treated as continuous. On such a finite volume, spatial momentum integrals 
are replaced by sums over the available momentum modes. Repeating the calculation of the loop 
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L [fm] l [fin] 



Figure 1: Volume dependence of the proton polarisabilities. Here we show the ratio of the difference of the 
finite and infinite volume results to the infinite volume results for three values of the pion mass using the 
parameters described in the text. The curves terminate at m x L = 3. 



diagrams using finite volume sums rather than integrals leads to the following expression for the 
loop contributions to 74: 

^° OP(L) = ~ U52np .L dX [ 3G B^(^)-4G T ^y 4 (^ A )\ , (3.4) 

where jtft = \Jm\ + X 2 and = y / m 2 +2XA + X 2 and the finite volume sums ^ a (m) and 
( m ) are defined in Ref. 

To illustrate the effects of the finite lattice extent, Figure [j] shows the volume dependence of the 
various polarisabilities in the proton. Here we have specialised to QCD, and chosen physical values 
for the various constants. In the figure, we show results for the ratio AX(L) = [X(L) — X(°°)]/X(oo) 
for the six polarisabilities at three different pion masses, m n = 0.25, 0.35, 0.50 GeV. The overall 
magnitude of these shifts varies considerably; generally volume effects are at the level of 5-10% 
for m n = 0.25 GeV and smaller for larger masses. Larger effects are seen in a number of the spin 
polarisabilities but there are as yet no lattice calculations of these quantities. Results in Ref. [|l|] also 
allow us to calculate the finite volume effects in the quenched data on the various polarisabilities 
calculated in Refs. [Q, Qj. The quenched expressions involve a number of undetermined LECs, so 
we can only estimate the volume effects. Using reasonable range for these parameters, we see that 
the calculations on a (2.4 fm) 3 lattice with m K =0.5 GeV may differ from their infinite volume 
values by 5-10% in the case of the electric polarisability. 
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4. Conclusion 

We have investigated Compton scattering from spin-half targets from the point of view of lat- 
tice QCD. We first discussed how external field methods can be used to probe all six polarisabilities 
of real Compton scattering for both charged and uncharged targets. We also discussed the effects 
of the finite volume used in lattice calculations on the polarisabilities. Since polarisabilities are 
infrared-sensitive observables (they scale as inverse powers of the pion mass near the chiral limit), 
they are expected to have strong volume dependence. This is indeed borne out in the explicit cal- 
culations presented here. In QCD, we generically find that the polarisabilities experience volume 
shifts of 5-10% from the infinite volume values for lattice volumes ~(2.4 fm) 3 and pions of mass 
0.25 GeV. In the case of quenched QCD (relevant to the only existing lattice data), we find signifi- 
cant effects even at pion masses ~ 0.5 GeV. Future lattice studies of the polarisabilities should take 
these effects into account in order to present physically relevant results. 
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